Atomistic Modeling of RuAl and (RuNi)Al alloys 
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Abstract 

Atomistic modeling of RuAl and RuAINi alloys, using the BFS method for alloys is performed. 
The lattice parameter and energy of formation of B2 RuAl as a function of stoichiometry and the 
lattice parameter of (Ru 5 Q. x Ni x )Al 5 Q alloys as a function of Ni concentration are computed. BFS- 
based Monte Carlo simulations indicate that compositions close to R^N^A^o are single phase 
with no obvious evidence of a miscibility gap and separation of the individual B2 phases. 


Introduction 


In comparison with nickel or cobalt aluminides, B2 RuAl has appreciable room-temperature 
toughness and plasticity and maintains considerable strength at high temperatures [1-2]. These 
properties, in combination with excellent oxidation resistance [3], make this alloy a potential 
candidate for aerospace applications though cost and high density are a significant concern. In an 
effort to drive down both cost and weight and improve upon its other properties, several studies 
[4-6] have looked at alloying schemes for replacing Ru or A1 with other elements that generally 
form an isostructural B2 phase such as Co and Fe for Ru and Ti for Al. But by far, the most widely 
studied ternary alloying addition has been Ni [7-13]. Even then, there is disagreement as to the 
structure of ternary Ni-Ru-Al alloys that exist between the NiAl and RuAl B2-phase fields. 
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Chakravorty and West [9-10] have reported a miscibility gap centered between the two binary 
phases resulting in a region consisting of two distinct B2 compounds. In apparent agreement, 
Sabariz and Taylor [11] have also observed a two-phase alloy at the composition Ni 2 5 Ru 2 5 Al 50 . 
While Homer et al. [13] found results similar to Chakravorty and West [10] in that man}, uui nui 
all, of the ternary compounds seemed to exhibit two distinct components, the evidence seemed to 
overwhelmingly suggest coring as opposed to actual formation of two distinct B2 phases. 
Furthermore, a sample within the miscibility gap claimed by Chakravorty and West was heavily 
milled and annealed so that diffusion distances would be much smaller and a better opportunity 
for obtaining a near equilibrium structure would exist. In this case, only a single B2 phase was 
observed. Furthermore, Liu et al. [12] observed only a single B2 phase across the (NiRu)Al 
system in mechanically alloyed samples when Ni was varied between 10 and 25 at.%, indicating 
complete mutual solubility between NiAl and RuAl. 

Given the lack of agreement for even the structure of Ni-Ru-Al alloys and the generally scarce 
data for other ternary or higher order systems, it is useful to consider the development of a 
modeling effort to supplement the ongoing experimental work on RuAl alloys. In this area, we 
provide the necessary tools to perform such modeling, and use them to study the fundamental 
properties of RuAl and (Ru,Ni)Al alloys. The modeling is performed within the framework of the 
BFS method for alloys, a quantum approximate method for the description of the energetics of 
complex systems at the atomic level, and is validated by comparison to experimental data. 

The BFS Method 

The BFS method for alloys [14] has been proven to be highly effective for the study of 
multicomponent systems. With the proper parameterization, it allows for an extremely 
economical, computationally simple, and physically sound description of the energetics of large 
collections of atoms. The BFS method is based on the assumption that the heat of formation, A H, 
of a given collection of atoms is the sum of the contributions of each atom in the sample, £j. Each 
contribution e, consists of two terms: a strain energy (e?) which accounts for the change in 
geometry with respect to a single monoatomic crystal of the reference atom, and a chemical 
energy (e, c ), linked by a coupling function (g,) so that £, = e, 5 + g,£, C . Three parameters, 
equilibrium lattice parameter, cohesive energy, and bulk modulus, for each of the constituent 


atoms are needed in the general derivative structure of the compound being studied (in this case, 
bcc) in order to calculate these terms. The chemical energy accounts for the corresponding change 
in composition, considered as a defect in an otherwise pure crystal. The chemical ‘defect’ deals 
with pure and mixed bonds, therefore, two additional perturbative parameters (A ab and A ba 
where A, B = Ru, Ni, Al) are needed to describe these interactions. All the parameters (listed in 
Table 1) are determined with the Linearized-Augmented Plane Wave method (LAPW) [15]. An 
alternative set of parameters, obtained with the Linear-Muffin Tin Orbital method [16] can be 
found in Ref. 17. We refer the reader to Ref. 14 for a detailed description of the BFS method and 
its parameterization, and Ref. 18 for details on Monte Carlo-Metropolis simulations using BFS. 


Results and Discussion 

A general comparison between the LAPW values and other calculated and experimental 
values of the lattice parameter, cohesive energy, and bulk modulus for RuAl [19-33] is shown in 
Table 2. The LAPW calculations correctly reproduce the value of the lattice parameter ag for 
stoichiometric RuAl. Although Ru and Al do not exist in the bcc phase, the lattice parameter of 
the computed bcc-Ru and bcc-Al are both larger than the equilibrium lattice parameter of RuAl (a 
1.4% difference with respect to bcc-Ru, and a 7.8% difference for bcc-Al). The predicted value of 
a 0 , 3.005 A, compares well with the experimental values, 2.95 A [32] to 3.03 A [20]. Considering 
the typical deviations in first-principles values for the energy of formation, the LAPW prediction 
compares well with the experimental value of 1.28 eV/atom [21]. The LAPW input, validated in 
Table 2, is used for the determination of the BFS perturbative parameters A ab and A ba , listed m 
Table 1 for all possible pair combinations of Ru, Al, and Ni. 

We first apply the parameters listed in Table 1 to the analysis of the B2 RuAl phase field, by 
considering a large number of computational cells (defined previously in Ref. 14) including all 
possible defect structures for off-stoichiometric compositions. Whether structural vacancies are 
included or not in this set of possible configurations, it is seen that the RuAl phase is quite 
unstable with respect to changes in stoichiometry, in the sense that there is a noticeable change in 
energy of formation off stoichiometry. Furthermore, this change is more pronounced for Al-rich 
alloys, as can be seen in Fig. 1. For Ru-rich alloys, the lowest energy configurations correspond to 
substitutional alloys. The analysis suggests a trend towards phase separation, as configurations 


showing clustering of excess Ru atoms are energetically favored. This can be explained in terms 
of the individual BFS contributions e, (i = Ru, Al) to the energy of formation of the computational 


cell. For B2 RuAl, e Ru = -0.3383 eV/ atom and e Al = -1 .6516 eV/atom, indicating that both atoms 


contribute to the formation of the alloy (negative contributions to AH). The ‘weak contribution of 


Ru atoms (compared to that of Al atoms) and the similarity between the atomic volume per atom 
in the B2 RuAl alloy and the atomic volume in a pure Ru crystal, lead to the phase separation 
observed for Ru-rich alloys [6]. To complete the analysis of B2 RuAl, Monte Carlo simulations 
using the NN approximation [18] were used to determine the melting temperature. It was found 
that the structure is highly ordered and stable up to 2300 K, in excellent agreement with the 
experimental value of 2330 K [1]. 

Given the possiblity of 1) reducing the cost and weight of RuAl alloys through substitution of 
Ru with Ni, 2) the potential to use RuAl as a reinforcing phase in other intermetallic systems, 
possibly even NiAl [7], and 3) given the controversy surrounding the phase structure of ternary 
Ru-Ni-Al alloys in the region between the two binary B2 phases NiAl and RuAl [9-13], Ru-Ni-Al 
would appear to be an ideal system for further theoretical analysis. While the following BFS 
analysis is applied to Ni additions to RuAl, it should be noted that similar calculations can be 
easily performed for other ternary or higher order additions using, for example, the BFS 
parameters listed in Ref. 17. 

For the ternary case, we first determine the preferred site occupancy for dilute additions of Ni 
in RuAl. Let A(B) denote an A atom occupying a site in the B sublattice, and A(B)B(C) the same, 
but with the displaced B atom moving to a site in the C sublattice. For (Ru 5 o- x Ni x )Al 50 alloys, 
Ni(Ru) = -0.98217 eV/atom and Ni(Al)Al(Ru) = -0.8780 eV/atom, indicating that Ni prefers 
available Ru sites. However, for Ru 50 (Al 50 . x Ni x ) alloys Ni(Al) =-0.9600 eV/atom and 
Ni(Ru)Ru(Al) = -0.9513 eV/atom, indicating a slight preference for Al sites, although the 
difference in energy is so small that it is difficult to conclude any obvious preference for either site 
under these conditions. The strong preference of Ni for Ru sites when Ni+Ru < 50 at.% is not 


surprising, as it allows for the formation of strong Ni-AI bonds. 

In terms of the phase structure of the ternary alloy(s) bridging the two binary B2 phases NiAl 
and RuAl, there are reports that suggest either a miscibility gap and phase separation for 
compositions near N^Ru^Aljo [9-11] or a single ternary B2 phase throughout [12-13]. While 
BFS Monte Carlo NN simulations, as described in Ref. 18, clearly show the formation of the B2 


RuAl phase, as can be seen in Fig. 2.a, the results are just as clear for N^Ra^AIso, as shown in 
Fig. 2.b. The simulations show that both Ru and Ni share the same sublattice. The degree and type 
of ordering can be estimated from the coordination matrices a and b, for nearest-neighbors (NN) 


and next-nearest-neighbors (NNN), respectively, where the matrix element a y 


n (h l denotes the 


probability that an atom of species i has an atom; as a NN (NNN). If there is phase separation, it 
would be expected that atoms in the Ru/Ni sublattice will have, mostly, atoms of their same 
species as NNN’s. If there was a ternary ordered phase other than B2, the NNN coordination of 
Ru or Ni atoms would reflect this fact with an increasing number of NNN’s of the other species. 


For example, the NNN coordination matrix for a Fleusler L2j ternary phase would maximize the 
value oib NiRu ( = b RuNi = 1). Phase separation, on the other hand, would lead to a very small value 
of b NiRu (and b RuNi ), maximizing in turn, b NiNi (and b RuRu ). Instead, the results of the simulations 
reflect the existence of a (Ru,Ni)Al B2 phase where Ru and Ni are randomly located throughout 
their own sublattice, in agreement with Liu et al. [12] and Homer et al. [13]. 

The BFS values of the lattice parameter for (Ru 50 _ x Ni x )Al 50 as a function of Ni concentration 
(in at%), a(x), show a slight deviation with respect to the the average values obtained using 
Vegard’s law [34], as shown in Fig. 3. Reported measurements of a(x) [1 1,12] are also considered 
as having a linear behavior. The linear fit of the experimental values in Ref. 1 1 (normalized to the 
equilibrium B2 value a 0 ) is a(x)/a 0 = 1 - 0.0006949x. In spite of the slight positive deviation from 
linearity, the best linear fit of the BFS prediction is almost identical, a(x)/a 0 = 1 - 0.0007266x. 

To further understanding of the behavior of the ternary alloys, it is possible to define within the 

framework of BFS, the role that each atom plays in (Ru 50 - x Ni x )Al 5 o as a function of Ni 

concentration, by separately computing the BFS energy contributions. Fig. 4 shows the strain, 

chemical and total energy contributions of Ru, Ni, and Al atoms to ternary RuNiAl alloys, 

c 

computed as the average over all atoms of similar species. Ru atoms provide a favorable (i.e., e < 
0, favoring alloying) chemical contribution to the total energy of formation, diminished by 
increasing strain as the Ni concentration increases. As a result, the role of Ru in favoring 
compound formation is diminished with increasing Ni content (Fig. 4.a). Ni atoms display the 
opposite behavior: as the concentration of Ni increases, the strain energy becomes substantially 
smaller as the average volume per atom in the alloy becomes closer to that of Ni. The chemical 
contribution, which does not favor alloying, increases, but at a slower rate than the decrease in 
strain, resulting in a total contribution that does not favor the stability of the alloy (Fig. 4.b). Al 



atoms display a mild increase in strain with increasing Ni content, easily compensated by the 
growing (much more favorable than Ru) chemical contribution, leading to a net decrease of the 
total energy (Fig. 4.c), which along with the chemical contribution from Ru ultimately favors the 
formation of ternary B^. ^Ouipvuiius. 


Summary 

The BFS method for alloys has been used to model the structure of RuAl and the extension of 
the B2 phase fields from NiAl to RuAl in the Ru-Ni-Al system. The BFS analysis correctly 
predicts the melting temperature of RuAl and indicates a strong preference for Ni in Ru sites in 
(Ru 5 o. x Ni x )Al 50 alloys, and an almost equal probability for either site in Al-defficient alloys. The 
results indicate that the two phases are isostructural with a single B2 phase of varying Ru:Ni ratio 
bridging the binary compounds. The calculated lattice parameter for (Ru 50 - x Ni x )Al 5 Q alloys, as a 
function of Ni contents, displays slight deviations from the average values, but in close agreement 
with experiment. 
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TABLE CAPTIONS 


Table 1: LAPW results for the lattice parameter, cohesive energy, and bulk modulus for the bcc 
phases of Ru, A1 and Ni, the resulting equivalent crystal theory (ECT) parameters p, a, / and X 
and the BFS perturbative parameters A ab and A ba . 


Table 2: Calculated (LAPW) properties compared to measured and other theoretical values for 
NiAl and RuAl. 


Table 3: Nearest-neighbor (NN) and next-nearest-neighbor (NNN) coordination matrices for the 
Ru-Ni-Al simulation shown in Fig. 2.b. The matrix element in row i and column j, denotes the 
probability that an atom of species i has a NN (top matrix) or a NNN (bottom matrix) of species), 
where i (or j) = 1, 2, 3 corresponds to Ru, Al, Ni, respectively. 


FIGURE CAPTIONS 


Fig. 1 : Energy of formation (in eV/atom) vs. Ru concentration. The data points shown correspond 
to the lowest energy configuration found at each composition. 


Fig. 2 : Final state (room temperature) of Monte Carlo MC-NN simulations for (a) RuAl and (b) 
Ni 2 5Ru25Al 5 o. Ru, Ni, and A 1 atoms are denoted with white, dark, and light grey spheres, respec- 
tively. 

Fig. 3 : Comparison of experimental results and BFS predictions for the lattice parameter of (RU50- 
x Ni x )Al 50 alloys as a function of Ni concentration. Experimental results are from Ref. 11 . The 
dashed line indicates the results of Vegard’s law. 

Fig. 4 : Individual average BFS contributions for a) Ru, b) Ni and c) A 1 atoms in a (Ru50-x^i x )Al50 
alloy, as a function of Ni concentration. In each case, the strain (Triangle up), chemical (triangle 
down) and total (solid disk) energies are shown. 


BFS parameters (in A' 1 ) 

Ar u ai -0.04186 

^RuNi 

0.28503 

A NiAl 

-0.04288 

AaiRu -0.02483 

A NiRu 

-0.04361 

A AlNi 

0.10287 


Table 1 




a 0 (A) 

B (GPa) 

A H (ev/atom) 

NiAl Expt. 

2.8864 [20] 

158-166 [23], 189 [24] 

0.61 [22], 0.6831 [25] 

LAPW 

2.89 

185 

0.79 

(This work) 

Other calcs. 

2.86 [26], 2.88 [27,28], 
2.839 [29] 

200 [26], 170[27], 
184 [29] 

0.697 [29], 0.68 [30], 
0.61 [19], 0.733 [31] 

RuAl Expt. 

3.03 [20], 2.99 [26], 
2.95 [32] 

208 [26] 

1.28 [21] 

LAPW 

3.005 

203 

0.99 

(This work) 

Other calcs. 

2.967 [29], 3.02 [31] 

223 [29], 220 [31] 

0.61 [29], 0.776 [31], 
0.81 [33] 


Table 2 


NN 


Ru 

A1 


Ru 

0.000 

0.992 


A1 

0.496 

0.016 


Ni 

0.008 

0.976 


NNN 

i/j 

Ru 

A1 


Ru 


0.008 


A1 

0.004 

0.985 


Ni 

0.512 

0.023 


Ni 

0.008 

0.488 

0.016 

Ni 

0.512 

0.011 

0.4653 


Table 3 







Fig. 2 








